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The four-dimensional helicity regularization scheme is often used in one-loop QCD computations. 
It was recently argued in Ref. [l[ that this scheme is inconsistent beyond the one-loop order in 
perturbation theory. In this paper, we clarify the reason for this inconsistency by studying the 
perturbative expansion of the vector current correlator in one-flavor QED through three-loop order. 
We develop a simple, practical way to fix the four-dimensional helicity scheme using the idea of 
dimensional reconstruction, and demonstrate its application in several illustrative examples. 



The four-dimensional helicity (FDH) regularization 
scheme 0, Q is one of several regularization schemes 0- 
6] based on the idea that consistent definitions of quan- 
tum field theories can be achieved through analytic con- 
tinuation in the number of space-time dimensions Q. 
The modern use of the FDH scheme is largely restricted 
to one- loop computations 0- The motivation for FDH 
arose from on-shell methods for loop computations Q, 
which seek to reconstruct higher-loop scattering ampli- 
tudes from tree-amplitudes through their unitarity cuts. 
The tree amplitudes in massless QCD have a remark- 
ably simple form Q if a four-dimensional concept - the 
spinor-helicity formalism fTp| - is employed in their eval- 
uation. If this simplification is to be used in loop compu- 
tations, the spin degrees of freedom for virtual particles 
must be treated as four-dimensional, in contrast to their 
momenta. This distinction is made manifest in the FDH 
scheme. 

Until recently, very little work was done to extend 
the FDH scheme beyond one-loop computations. To the 
best of our knowledge, the majority of complete multi- 
loop computations in non-supersymmetric theories are 
performed using conventional dimensional regularization 
(CDR) (2t|. In contrast to this, the known higher-order 
FDH results include some two-loop scattering amplitudes 
(see Ref. [3| for examples), that have not been used for a 
computation of any physical quantity. 
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There is a good reason for this state of affairs. The 
CDR scheme is a natural scheme to use if inclusive quan- 
tities such as cross-sections and decay rates are computed 
using the optical theorem. Since computations of multi- 
loop integrals in those cases are mostly based on the 
integration- by-parts identities (Til Il2j | , it is important to 
set up calculations in such a way that D-dimensional 
Lorentz invariance, where D is the space-time dimension 
of CDR, is explicit in all stages of the computation. The 
fact that intermediate observable states are more natu- 
rally described by four-, rather than D-dimensional quan- 
tum fields, is immaterial within such an approach. The 
distinction between observable and unobservable states 
is accomplished indirectly, by taking the imaginary part 
of an appropriate Green's function at the very end of the 
calculation. 

However, CDR is also used in existing fully differential 
next-to-next-to-leading order (NNLO) computations [T3T - 
[l9l ] . Its use in such situations is much less natural. CDR 
necessitates calculations of multi-parton matrix elements 
to higher orders in e = (4 — D)/2. If FDH were extended 
to the two-loop order, this step could be avoided, leading 
to increased efficiency in computations of quantities with 
rich multi-parton kinematics. In fact, FDH is a scheme of 
choice in many calculations that address next-to-leading 
order (NLO) QCD corrections to kinematic distributions 
at hadron colliders (see Ref. 20] for a recent summary of 
results). Having in mind that extension of perturbative 
computations for some basic LHC processes to NNLO is 
desirable for several reasons, it is interesting to under- 
stand if the FDH scheme can be used beyond one loop 
in non-supersymmetric theories. 
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In a recent work [l[ , Kilgore made a step in this direc- 
tion by studying the application of several regularization 
schemes to higher-order calculations: CDR, FDH, and 
the dimensional reduction approach commonly used 
in supersymmetric theories. He considered the imaginary 
part of the correlator of two vector currents and pointed 
out that the FDH scheme as formulated in Ref. Q be- 
comes inconsistent at higher orders. 

The goal of this article is to elucidate the reasons be- 
hind this inconsistency, and see if they can be fixed. We 
find that the problem with the FDH scheme follows from 
the fact that the gauge invariance of the full theory is bro- 
ken by the restriction of the loop momenta to a smaller 
dimensionality than the spin dimensions. As the result 
of this, the Ward identities are not satisfied for the ad- 
ditional spin degrees of freedom. The situation with the 
FDH scheme becomes very similar to that of the dimen- 
sional reduction. However, we also find that the ills of 
FDH at NNLO can be cured very simply using "dimen- 
sional reconstruction" : if the one-loop result for an ob- 
servable is known for an arbitrary number of spin di- 
mensions, then the (incorrect) two-loop FDH result can 
be fixed once the one-loop renormalization constants of 
the theory are known for two different integer numbers 
of spin dimensions. No two-loop CDR computation is 
necessary, so that this set-up preserves some of the sim- 
plicity of the original FDH. The spirit of this fix is similar 
to a technique employed to reconstruct the rational parts 
of one-loop amplitudes in generalized D-dimensional uni- 
tarity plj . 

We would like to illustrate this idea by considering as 
simple a set-up as possible. We choose to study Quantum 
Electrodynamics (QED) with a single massless fermion 
field. The Lagrangian of the theory reads 



L = i>{id^ + eA^)*ij> 



4 ^ ' 



(1) 



where F^ v = d^A v — d v A^. The vector current in this 
theory J M = ■tpj^ip is conserved, d^J^ = 0, and does not 
require renormalization. We study the correlator of two 
vector currents, 



^(q)n(q 2 ) = -il d 4 xe^ x (0\TJ^x)M0)\0), (2) 

where we used uj^,j(q) = {q 2 g^i, — q^q v )- It is important 
to stress that, as it is customary in FDH, the external 
Lorentz indices fi and v are taken to be four-dimensional. 
Conservation of the vector current requires that the imag- 
inary part of the correlator Im[n(g 2 )] is finite without 
any renormalization if the gauge-invariant subset of di- 
agrams with closed fermion loops (the only contribution 
to the renormalization of the electric charge e) are dis- 
carded. It is this absence of any renormalization that 
makes this quantity an ideal laboratory for our investi- 
gation into various regularization schemes. We note that 
Kilgore studied the correlator of the two conserved vec- 
tor currents in QCD [l| , where disentangling the coupling 
constant renormalization is possible, but more difficult. 



We believe that focusing on the QED aspect of the prob- 
lem allows us to illustrate the main issue very sharply. 

It is straightforward to compute H(q 2 ) to three- loops 
using various regularization schemes. We have done this 
in the variety of ways, including utilizing Mincer 22] or 
Air [23j], as well as using in- house implementations of 
the Laporta algorithm [2J| established previously [25j . 
We find complete agreement with the result reported 
in Ref. |l| when it is truncated to QED and all terms 
that are proportional to the number of lepton flavors are 
dropped. The imaginary parts of the correlator com- 
puted in CDR and FDH are 
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The difference is striking. It implies that the computa- 
tion of a finite quantity, that does not require any renor- 
malization, leads to different results when two different 
regularization schemes are applied. Moreover, we empha- 
size that both CDR and FDH computations are consis- 
tent with the conservation of the vector current J M in four 
dimensions so there is nothing at this point that makes 
either of the two results in Eq. (J3)) obviously incorrect. 
One could have suspected that a finite shift in the cou- 
pling constant - familiar from the application of the FDH 
scheme in one-loop QCD computations Q- can account 
for the difference of the two results. However, the known 
shift of the coupling constant is purely non-abelian Q- 
Since it vanishes in the abelian (QED) limit, it is not 
possible to reconcile the two results shown in Eq. ^ by 
existing means. 

To understand the reason behind the difference, we re- 
view the rules Q that are used in the FDH computation. 
We begin by considering QED in a D s -dimensional space, 
but with all momenta restricted to a /^-dimensional sub- 
space of this Z? s -dimensional space. This arrangement 
requires D s > D. Upon performing spin algebra in all 
contributing diagrams, we take D s — » 4, keeping D fixed. 
The limit D — > 4 is taken at the end of the calculation. 
We note that the CDR scheme can be formulated in a 
similar way, making it explicit that the two schemes dif- 
fer by the order of limit-taking. Indeed, to arrive at the 
CDR result, we take D s D for fixed D, and then take 
the limit fl44. 

The origin of the differences in CDR and FDH results 
can be best understood by presenting Tl(q 2 ) in a form 
where D s is kept fixed, while the limit D —> 4 — 2e is 
taken. We find 
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where S s = Dg—A. The FDH result is obtained by setting 
6 S = in Eq. while the CDR result corresponds to 
setting S s = — 2e and taking the limit e — > 0. 0(S S ) terms 
that are not enhanced by inverse powers of e are present 
at NNLO but are irrelevant for both CDR and FDH. A 
similar term at NLO is also irrelevant for both CDR and 
FDH but, as we will see, it is important for understanding 
differences at NNLO between the two schemes. For this 
reason, it is shown explicitly in Eq. (|4|. 

It follows from Eq. Q that the difference between the 
CDR and FDH schemes appears at NNLO because terms 
of the form a 2 S s /e are present in that order of the per- 
turbative expansion. Those terms either contribute to 
the final result (CDR), or are set to zero by convention 
(FDH). Note that no 8 s /e term appears at NLO. There- 
fore, to understand the difference between CDR and FDH 
schemes, we must explain why divergent terms propor- 
tional to the number of "extra-dimensional" degrees of 
freedom appear at NNLO. 

The reason becomes very clear if we set D s to an in- 
teger value greater than four. For the sake of argument, 
we take D s = 5. It immediately follows from Eq. (U) that 
ImII(g 2 ) is divergent. To see why this divergence occurs, 
we must go back to the QED Lagrangian in Eq. ((T|) and 
ask what theory arises if we set D s to five but keep all 
space-time coordinates four-dimensional. 

We begin by extending the Dirac algebra to five di- 
mensions by taking r M = 7 M , fi = 0, .., 3 and T 4 = 275, so 
that 



pJVpM _ 2g MN 



M,N = 0,.A. (5) 



The fermion fields are not analytically continued, so the 
number of independent fermion helicities remains two. 
The gauge field A is split into a four-dimensional gauge 
field and a scalar field, A M = {A»,4>). The QED La- 
grangian of Eq. (TTJ) written in terms of four-dimensional 
fields reads 



(G) 



Note that in Eq. © we introduced a new coupling con- 
stant <?0, to parameterize the interaction of the field <f> 
and the pseudoscalar fermion current. Because the La- 
grangian in Eq. © originates from the five-dimensional 
QED Lagrangian in Eq. (TTJ), = e. However, since 
we use four-dimensional momenta and coordinates, this 
equality of the coupling constants can not be protected 
by the full D s -dimensional gauge invariance. This im- 
plies that in Eq. (|6]), the coupling constant g^ requires 
renormalization, while the electric charge e is not renor- 
malized and is protected by the four-dimensional gauge 
invariance. For D s = 5, the result shown in Eq. ^ 
corresponds to the calculation of Im[n(g 2 )] in a theory 
defined by Eq. © in terms of bare charges, e and g^. 
The bare electric charge e coincides with the physical 
charge because of the four-dimensional gauge invariance, 



but renormalization is required for g^ to make the corre- 
lator of the two vector currents explicitly finite in D s = 5. 
Such renormalization has not been performed in Eq. (Q| . 
This is the reason for the 1/e divergences present there. 
We conclude that the "divergences" in Eq. Q - crucial 
for understanding the CDR/FDH difference - can be re- 
lated to the renormalization of the coupling constant g^ 
for finite D s . Below we describe the details of this rela- 
tion. 

Because the scalar field 4> contributes to the corre- 
lator of the two vector currents only at NLO, through 
NNLO we only need the one-loop renormalization of the 
coupling constant g^. It is easy to obtain this renor- 
malization constant by considering the Green's function 
(0|TV;(a;)^(x)|0). We find 
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where = g^f^n is introduced. Rewriting Eq. ((4]) for 
D s = 5 (S s — > 1) and separating the two couplings at 
NLO explicitly, we find 



Im[n( g 2 )]^ =5 = I i- 
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Removing the bare coupling from Eq. © using Eq. (J7J, 
we see that the divergence in Eq. ([5} disappears. This 
proves our assertion about the origin of the divergent S s /e 
terms in Eq. (@]). 

Having understood the origin of divergences in Eq. ((4]) , 
we must find a way to calculate the difference between 
Im[n(g 2 )] in the FDH and CDR schemes without per- 
forming a complete three-loop computation. We observe 
in Eq. (|4]) that only the 0(S s /e) term contributes to the 
CDR/FDH difference; the 0(<5 2 /e) term is not relevant. 
However, since Eq. (HJ) is a second-degree polynomial in 
S s , it is not possible to isolate the desired term by per- 
forming the computation in a single integer-dimensional 
space. Two such calculations are required. A similar 
need occurs when attempting to reconstruct the rational 
parts of one-loop amplitudes using tree-level amplitudes 
in higher integer dimensions [2l| . albeit for a different 
reason. 

We have already discussed the case D s = 5. The case 
D s = 6 is qualitatively similar, but different in detail. 
The D s — 6 QED Lagrangian of Eq. |T]) deconstructs to 

L =^(i llM d^ + eA M 7^)* - g4> V2 (^ 75 ct+*0 + h.c) 

1 (9) 

where ^ = (u, d) is the lepton "doublet" , a + — (a 1 + 
ia 2 )/2, a 1,2 ' 3 are the Pauli matrices and is a complex 
scalar field. We define the conserved vector current as 
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J/j, = (ujnU + df) fj,d) I 'V2) where the normalization factor 
is chosen for convenience. The corresponding result for 
the imaginary part of the polarization operator follows 
from Eq. where we isolate the contribution due to 
scalar degrees of freedom at one-loop: 



im[n(g 2 )]„ =6 



1 
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The divergence is removed by the renormalization of the 
bare coupling g^ which, for D s — 6, is computed from 
the "flavor-changing" Green's function (0\Td(f)u\0) . We 
find 



, bare I f i 
*U \D S =6 = OC I 1 
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It is clear from Eq. fTU]) that this renormalization of the 
coupling constant makes Im [n(g 2 )] finite. 

Since we understand the structure of ultraviolet diver- 
gences for two values of D s , it is easy to find a rela- 
tion between the FDH and CDR schemes. We imagine 
that a one-loop computation is performed, and the In- 
dependence of the one-loop result is established. We as- 
sume that the two-loop FDH result is also known. The 
result for general D s reads 



im[n( g 2 )] 5s =im[n( g 2 )] FDH - 
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The CDR result corresponds to setting S s — > — 2e in 
Eq. (p~2j) and neglecting all 0(e) terms. Doing so, we 
find 



Im[n(a 2 )] CDR = Im[n(g 2 )] 



fdh cwoV (13) 



The connection between the two schemes requires knowl- 
edge of the coefficient Ci. As we discussed earlier, both 
ci and C2 are related to the renormalization constants of 
the couplings of pseudoscalar fields, that appear as the re- 
sult of dimensional deconstruction, to fermion bi-linears. 
Hence, it is a simple matter to find c\. We require that 
Eq. (TT2"T) becomes finite for D s = 5, 6 if the renormaliza- 
tion of the coupling constant is performed. Doing so 
for both values of D s leads to a system of two equations 
that can be solved for c\ and C2- We find 
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In Eq. (|14[) . we have introduced the one- loop renormal- 
ization constants for the couplings of the scalar fields to 
fcrmions in the compactification of Z3 s -dimensional QED 
to four-dimensional space-time: 



^ bare I 
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(15) 



Explicit expressions for SZ[j s for D s = 5,6 follow from 
Eqs. (|7llip . Using those results, we find 

12rfm [n(g 2 ) CDR - H(o 2 ) FDH ] = | ) 2 , (16) 

in agreement with the explicit computations of Eq. ([3]) . 
As advertised, we are able to obtain the correct NNLO 
result from the FDH result without dealing with D s - 
dimensional spin degrees of freedom at NNLO. 

There are several possible directions that one can ex- 
plore at this point, including how this picture generalizes 
to more complicated theories (QCD, massive QED, etc.) 
or more complicated observables. Except for a few com- 
ments, in this paper we restrict ourselves to QED but we 
study observables that depend on the mass of the lep- 
ton. We show that the procedure we introduced in the 
context of the vector current correlator is general and re- 
mains valid also in those cases. We work with one massive 
fermion flavor in both examples. 

• We begin by computing the mass renormalization con- 
stant in FDH in the on-shell scheme, and ask if we can 
relate it to the mass renormalization constant in CDR. 
The mass renormalization constant is defined as 



m 



(17) 



where mo is the bare fermion mass and m is pole mass of 
a lepton. One can easily read off the MS mass renormal- 
ization constant from Eq. (fTTf because the lepton pole 
mass is an infra-red finite quantity. We compute Z m 
through two-loop order in QED. We consistently neglect 
the contribution of the fermion loops, so that the elec- 
tric charge does not need to be renormalized. The mass 
renormalization constant takes the form 

Z m = l + tt „zW+o^>, (18) 

where ao = a/7rr(l + e)/(47r)~ e TO~ 2e and 
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4e 4 8 V e 



+ 1 +0(e), 



m e 2 ^32 + 16 128/ t ^64 16 1 J 

13S 2 S \ 219 5tt 2 tt 2 , n 3 A m , . 

1 H 1 In 2 ( 3 + 0(e). 

256 J 128 16 2 4 Si v ; 

We have written the result in a form where the (In- 
dependent terms are manifest. The expression Eq.(|19p 
can be translated into the CDR and FDH values for 
the on-shell mass renormalization constant by setting 5 S 
to the appropriate values. Suppose we have computed 
Eq. (TT9")) using the FDH scheme. Can we obtain the mass 
anomalous dimension in the CDR scheme without doing 
a complete calculation? 

The evolution equation for the mass parameter reads 



dm(fi) 



m (2ea 
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Taking Z™ H from Eq. ([19]) and setting /3(a) 
find 



d/i 



7717(a) 



which implies 
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To find the anomalous dimension in the CDR scheme, 
we write a relation between the FDH renormalization 
constant and the renormalization constant at arbitrary 



Z m (8 s )-Z, 
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where the ellipses stands for other terms that do not af- 
fect the anomalous dimension. To find Z^ DR , we need 
c 2 i, since it leads to divergent contribution in the limit 
8 S = — 2e. Repeating what we did for the photon vacuum 
polarization, we must consider the theory at finite D s , so 
that C21 and C22 contribute to the leading two-loop diver- 
gence of the fcrmion self-energy Since such divergence is 
entirely fixed by the lowest-order mass anomalous dimen- 
sion and the /3-functions for the coupling constants, we 
can find an equation for c 2 i and c 22 . We note that the 
/3-functions appear because of the need to renormalize 
the scalar-fermion couplings, as described in Eqs. (|7|11|) . 
The relevant condition is that 



d d 
2ea— \nZ m (6 s ) + ,S(cu)- — \nZ m (5 s 
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(24) 



is free from 1/e singularities for any value of S s . In prac- 
tice, we choose D s — 5 and D s = 6. The /3-functions 
follow from Eqs. (|7|lip . We write them here for com- 
pleteness: j3{a$) = — 3/4a 2 for D s = 5, and /3(a$) = —a 2 
for D s = 6. We finally find c 2 i = 1/16 and c 22 = -1/128, 
in agreement with Eq. (|19j) . The mass anomalous dimen- 
sions in the CDR scheme follows immediately. Finally, 
one can imagine that the difference between on-shell Z m 
factors in different schemes Eq. ([T9l) can be understood 
completely, by going beyond the MS renormalization of 
the g$ coupling constants as in Eqs. (171111) and insisting 
that the two couplings g$ and e are equal to each other, 
including the finite renormalization. We did not pursue 
this question in this paper but it is an interesting avenue 
for further studies. 

• As the final example we compute the two-loop QED 
corrections to the electron anomalous magnetic moment 
and show that the correct result can be obtained using 
the FDH scheme and the procedure outlined above. We 
begin by writing the amplitude for the electron scattering 
off the electromagnetic field as 



iM = —ieu(p2)Tu(pi), 
T = eFAQ 2 ) + ^Z^F2(q 



(25) 



In Eq. (|25[) , e M is the "polarization vector" of the external 
field, e = 7^6^, and q — p 2 — pi is the momentum transfer 
from the electron to the field. The anomalous magnetic 
moment is given by a e — (g — 2)/2 = F 2 (0). The one-loop 
result for arbitrary 5 S is given by 



(26) 



The two-loop result for g — 2 requires the on-shell wave- 
function and mass renormalization constants for the elec- 
tron at the one-loop order. The mass renormalization 
constant Z m is given in Eq. ()19[) . The wave-function 
renormalization constant Z 2 coincides with Z m in QED 
at this order in both CDR and FDH schemes. We find 
the following results for the two-loop contribution to the 
electron anomalous magnetic moment in the CDR and 
FDH schemes 
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Our CDR result matches well-known results in the liter- 
ature [26l . [27| , when fermion-loop contributions are ne- 
glected. The FDH result is new. We now illustrate how 
to use dimensional reconstruction to obtain the CDR re- 
sult, given the (^-dependent 1-loop result in Eq. (|26| and 
the 2-loop FDH result. We proceed as we did for the cur- 
rent correlator by writing the result for arbitrary S s as 
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The CDR result is obtained by taking 5 S 



-2c. 
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To obtain ci, we compute Eq. ([28)1 for D s =5,6 
demand the result be finite after renormalizing a^ lre . 
obtain 



(29) 

and 
We 
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Inserting this into Eq. ([29]) . we derive the correct (CDR) 
result for g—2. Hence, the procedure that we developed 
by studying the correlator of two conserved currents ap- 
pears to be valid in a more general context. 

Before concluding, we comment on two possible venues 
for the extension of this analysis, namely its extension 
to QCD and to its application to less inclusive observ- 
ables. The first comment concerns the well-established 
procedure for applying the FDH scheme in one-loop QCD 
computations. As explained in Ref. [7], it is possible to 
use FDH in one-loop computations consistently provided 
that a finite renormalization of the strong coupling con- 
stant, 
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is performed. We can easily understand this result us- 
ing dimensional reconstruction idea. Dimensional recon- 
struction in QCD leads to the appearance of color-octet 
massless scalars that interact with both fermions and 
"four-dimensional" gluons. Tree-level computations in- 
volve four-dimensional fields by definition, and therefore 
all one-loop amplitudes are proportional to the "four- 
dimensional" version of the strong coupling constant 
a s . Massless QCD is made finite by the coupling con- 
stant renormalization which, in the case of dimensional 
reconstruction, involves the contribution of color-octet 
scalars. Because we only need the divergent contribu- 
tion of massless color-octet scalar to the renormalization 
of a s , we can find it by inspecting the QCD (3- function, 
j3o = 11/3Ca — 2/3-/V/ — Ca/6N s and focusing on the con- 
tribution of the color-octet scalars (the term proportional 
to N s ). As expected, the required shift in the coupling 
constant in Eq. (j3Tj) and the contribution of the color- 
octet scalars to QCD /3-functions are appropriately cor- 
related. By studying the FDH scheme in one-loop QCD 
in terms of dimensional reconstruction, it is obvious that 
finite renormalization of the coupling constant in Eq. pip 
is the only thing needed to perform self-consistent com- 
putations in FDH 30]. 

As a second comment, it is interesting to ask what 
the dimensional reconstruction procedure outlined in this 
paper implies for exclusive computations. For the sake 
of argument, consider again the correlator of two vector 
currents in QED. Its imaginary part is directly related 
to the inclusive decay rate of a vector boson. But how 
should a decay rate be treated if we require a certain 
number of "jets" , borrowing from the QCD terminology? 
At NNLO, it is possible to have four, three and two jets 
in the final state. The four-jet rate is finite at this or- 
der. The three-jet rate is only needed through NLO, 
and therefore FDH can be used straightforwardly. The 
two-jet rate is needed at NNLO, which makes it obvious 
that this is the place where corrections to the inclusive 
rate must be accommodated. Moreover, since the phase- 
space for the two-jet configuration can be driven arbi- 
trarily close to the two-parton kinematics by appropriate 
adjustments in the jet selection criteria, the correction to 
the inclusive cross-section is the finite renormalization of 
the Born matrix element. This argument applies to pro- 
cesses which possess infra-red finite total cross-sections, 
but it needs further refinements for consistent applica- 
tion of the FDH scheme to exclusive processes at hadron 
colliders beyond one-loop. 



To conclude, in this paper we explored the four- 
dimensional helicity scheme at NNLO, following an in- 
teresting observation in Ref. [l[ that it becomes inconsis- 
tent at that order in perturbation theory. To avoid the 
complications of renormalization, we studied QED cor- 
rections to the imaginary part of the correlator of two 
conserved currents. We found that the differences be- 
tween the FDH and the CDR schemes are related to the 
fact that, upon continuing QED to a space-time of higher 
dimensionality while restricting all the loop momenta to 
lower-dimensional space-times, (D s — 4) components of 
the gauge fields turn into scalar fields and become un- 
protected by full D s -dimensional gauge invariance. As a 
result, divergences are introduced that require additional 
renormalization. They are removed in the FDH scheme 
by simply ignoring these additional degrees of freedom. 
In the CDR scheme, terms of the form (D s — 4) /e give 
additional finite contributions. One can argue for the 
correctness of the CDR result over FDH result by sim- 
ply stating that the former respects gauge invariance of 
the theory in D s -dimensions, while the latter only re- 
spects four-dimensional gauge invariance. Restoring the 
_D s -dimensional gauge invariance in the FDH scheme is 
possible using dimensional reconstruction: if the one- 
loop result is known for arbitrary S Sl and the two- loop 
FDH result is known, then the two-loop result in the 
CDR scheme can be obtained by studying the one-loop 
renormalization of the theory in D s = 5,6. This gives 
a simple, practical prescription for maintaining the sim- 
plifying features of FDH while still getting the answer 
right [3l[ . We demonstrated the applicability of this pro- 
cedure by reconstructing the two-loop CDR results from 
FDH for three examples: the correlator of two vector 
currents in QED, the mass anomalous dimension of a 
fcrmion in QED, and the electron anomalous magnetic 
moment. Further studies are required to extend these 
ideas to QCD and develop them to a point where com- 
putations of fully exclusive observables through NNLO 
in the FDH scheme are possible and transparent. This 
remains an interesting problem for the future. 
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